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ABSTRACT 

We present time series photometry for six partial transits of GJ436b obtained with the Fine Guidance Sensor instrument on the 
Hubble Space Telescope (HST). Our analysis of these data yields independent estimates of the host star's radius R+ = 0.505*2"^ Rq, 

and the planet's orbital period P = 2.643882^}°™!?° d, orbital inclination i = 85.80° + ° 2 ^, mean central transit time T c = 

2454455.279241+°-™^ HJD, and radius R p = 4.W°f 3 R B . The radius we determine for the planet is larger than the previous findings 
from analyses of an infrared light curve obtained with the Spitzer Space Telescope. Although this discrepancy has a 92% formal sig- 
nificance (1.7cr), it might be indicative of systematic errors that still influence the analyses of even the highest-precision transit light 
curves. Comparisons of all the measured radii to theoretical models suggest that GJ436b has a H/He envelope of ~ 10% by mass. We 
point out the similarities in structure between this planet and Uranus and Neptune and discuss possible parallels between these planets' 
formation environments and dynamical evolution. We also find that the transit times for GJ436b are constant to within 10 s over the 
1 1 planetary orbits that the HST data span. However, the ensemble of published values exhibits a long-term drift and our mean transit 
time is 128 s later than that expected from the Spitzer ephemeris. The sparseness of the currently available data hinders distinguishing 
between an error in the orbital period or perturbations arising from an additional object in the system as the cause of the apparent 
trend. Assuming the drift is due to an error in the orbital period we obtain an improved estimate for it of P = 2.643904 ± 0.000005 d. 
This value and our measured transit times will serve as important benchmarks in future studies of the GJ 436 system. 

Key words, techniques: photometric - eclipses - stars: individual: GJ 436 - planetary systems 



1. Introduction 

Transiting exoplanets present a unique opportunity to evaluate 
planet formation and evolution models. The principle informa- 
tion obtainable from the study of such planets are their masses, 
radii, and orbits. This complete data set cannot be obtained 
for non-transiting planets using any other methods with current 
technology. The radius of a planet for a given mass and inso- 
lation is a strong constraint for theoretical models and yields 
information about the planet's structure. Additionally, carefully 
designed observations of transiting exoplanets can yield data that 
were, up until recently, widely thought to only be the purview of 
future direct imaging studies. Such data include, but are not lim- 
ited to, exoplanets' chemical compositions, thermal emission, 
albedos, and energy redistributio ns (for a recent review of the 
field see Charbonneau et al. 2 007b . 

The recent discovery bv lGillon et alj (l2007bl) that the planet 
GJ436b, w hich was or i ginal ly identified with Doppler spec- 
troscopy bv lButler et al.l d2004l) . transits its host star is an impor- 
tant advance in this rapidly burgeoning field. It is the first mem- 
ber of the emerging class of planets known as "Hot Neptunes" 
(i.e. planets having masses similar to that of Neptune and or- 
bital semimajor axes ~ 0.01 - 0.10 AU) that has been observed 



* Table 2 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



to transit. Therefore, detailed observational study of an ice giant 
outside our solar system can now be carried out. 

Since the reported detection of transits for GJ436b there 
have been a flurry of follow-up observations and analyses of 
the system. So far, published results include complete analyses 
(i.e. data reduction and light curve modeling) by two separate 
groups of a single transit light curve and secondary eclipse light 
curve obtained at 8 pm using the InfraRed Array Camera (IRAC) 
on the Spitzer Spac e Telescope (hereaft er Spitzer for short). 
iGillon et all (l2007al hereafter G07) and iDeming et ail d2007l 
hereafter D07) analyzed the transit light curve and determined 
a precise mass, radius, and orbit for the planet. Both groups 
obtained consistent results and both concluded that the planet 
must have a gaseous H/He envelope similar to Neptune based on 
comparisons of its radius with theoretical models. Additionally, 
Southworth| (120081 hereafter S08) has recently analyzed the re- 
duced Spitzer light curve presented by G07 to give yet another 
estimate of the planet's radius, but did not use the result to make 
deductions about the planet's structure. The radius of GJ436b 
given by S08 is smaller than that ascertained by G07 and D07, 
although the results are consistent within the quoted uncertain- 
ties. 

D07 and lDemorv et al.1 d2007l) carried out complete analyses 
of the Spitzer secondary eclipse light curve for GJ 436b and were 
able to measure the planet-to-star flux ratio, and thus deduce the 
brightness temperature of the planet at the bandpass of the ob- 
servations. The two groups found nearly identical results, but 
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differed in their interpretation. D07 suggested that the measured 
temperature was higher than expected from a simple model of 
the planet and was indicative of tidal heating arising from the 
planet's non-circular orbit. Comparing the measured brightness 
temperature to more sophisticated models, Dem orv et alj 02007) 
found that a model with inefficient redistribution of energy from 
the day to night side could reproduce the observations. 

In this paper we present an analysis of a transit light curve for 
GJ436b obtained with the Fine Guidance Sensor (FGS) instru- 
ment on the Hubble Space Telescope (HST). Following the an- 
nouncement of transits for GJ436b we were granted Director's 
Discretionary time to carry out these observations. The motiva- 
tion was t o use the exception al photometric capabilities of the 
HST (e.g. lBenedict et al.ll 19981) to make a precise, independent 
determination of the radius and orbital ephemeris for the planet. 
Such measurements are of critical importance because they are 
the foundation on which the illuminating follow-up studies that 
are only possible for transiting planets are based. In the case 
of GJ 436b specifically, an independent measurement of the ra- 
dius is a test of previous assertions that the planet has a gaseous 
H/He envelope like Neptune and allows robust interpretation 
of secondary eclipse measurements. Additionally, the measured 
ephemeris will be a benchmark for future transit timing studies 
aimed at detecting additional planets in the system. 



2. Observations 

Six observations of partial transits for GJ 436b were carried out 
between 4 December 2007 and 5 January 2008 using the FGSlr 
on HST. This period spanned 1 1 orbits of GJ 436b. A log of the 
observations is given in TableQ] The FGSs are white light shear- 
ing interferometers that use four photomultiplier tubes (PMTs), 
two in each orthogonal direction on the sky, to precisely mea- 
sure tilts in the incoming wavefront. The measured counts from 
the PMTs are continuously recorded at a frequency of 40 Hz 
(25 msec). The instrument can be used as a high-speed photome- 
ter because of this, in addition to its more recognized uses for 
high precision relative astrometry and high angular resolution 
sensing. 

The FGSlr was used in "Position" mode for our observa- 
tions. In this mode the instrument uses the measurements of 
the wavefront tilt to lock on to the observed photocenter and 
hold fixed pointing. At a given time the FGS measures the 
photometric signal in a selectable 5"x 5" sub-field out of the 
greater instrument field of view (FOV). Position mode observa- 
tions do not allow simultaneous measurements of multiple tar- 
gets. Therefore, observations of multiple targets must instead be 
carried out for each one separately. 

We used the F583W filter for all observations. This filter has 
a nearly constant 95% transmission function over the range 450 
- 700 nm, and drops sharply to 0% outside that range. The sensi- 
tivity of the PMTs decreases linearly from a maximum of 15% to 
4% towards redder wavelengths over the non-zero transmission 
range of the filter. 

The observations consisted of six visits with a single orbit 
each. During each orbit we observed two reference stars with 
similar brightnesses (GSC 1 984-02050 and GSC 1 984-02056) 
as GJ 436 before and after the main observation of it to facilitate 
correcting for inter- and intra-orbit variances in the photome- 
try. The length of each reference star observation was 100 s. The 
main observation of GJ436 in each orbit was 1800 s. The ob- 
serving sequence for all visits was: reference star #1, reference 
star #2, GJ 436, reference star #2, and then reference star #1 . 



Table 1. Log of GJ436 observations. 



UT date Orbit phase (T c ) Relative flux level 



4 Dec 2007 
7 Dec 2007 
12 Dec 2007 
28 Dec 2007 
3 Jan 2008 
5 Jan 2008 



0.9942 - 0.0021 
0.0016-0.0095 
0.9917-0.9996 
0.9866 - 0.9945 
0.0015 -0.0094 
0.9838-0.9918 



0.99962 ± 0.00023 
0.99967 ± 0.00017 
0.99906 ± 0.00021 
1.00358 ±0.00011 
0.99782 ± 0.00021 
1.00024 + 0.00010 



Initial processing of the raw data files was done with the FGS 
data reduction pipeline. This yielded the measured counts from 
each of the four PMTs, the calculated position of the target in the 
FOV, and data quality flags for each individual 40 Hz sample. We 
retained readings obtained only when the instrument was in fine 
lock on one of the targets and no errors were indicated. 

The pipeline also output the average counts per 40 Hz sample 
for each of the PMTs when moving the instantaneous FOV from 
one target to the next during an orbit. This is a good estimate 
of the total dark and background counts for the PMTs. The esti- 
mated values for all four PMTs during all si x visits agreed with 
the res ults from the nominal values given by iNelan & Makidorj 
(120021) . The values were on average ~ 0.2% of the measured 
counts for GJ 436. The first step of post-pipeline processing was 
then to subtract the corresponding background from each 40 Hz 
sample measured by the four PMTs. 

After the background subtraction, we corrected each sample 
for the dead time i n the PMTs according to the prescription of 
INelan & Makidon] 12002). We converted the geocentric observa- 
tion times to the heliocentric frame using the recorded position 
of the spacecraft relative to the Sun and the observation vector 
to GJ436. The correction from the spacecraft frame to the geo- 
centric frame was negligible. 

During observations in fine lock on a target the FGS contin- 
uously corrects the fine pointing using the measured count dif- 
ferences in the two orthogonally oriented PMT pairs. In an ideal 
instrument, this difference would be zero when the instrument 
was pointed directly at a target. However, the PMTs have dif- 
ferent sensitivities and a parallel wavefront will result in small, 
but important, differences in PMT pair counts. This wouldn't be 
a problem for photometry in and of itself, except that the fine 
pointing of the instrument is not perfect and the apparent posi- 
tion of the star varies randomly. The impact of these two issues 
is that neither the average nor the total of the PMT counts in each 
pair is a conserved quantity. To correct for this we estimated the 
relative sensitivities of the PMTs in each pair using the aver- 
age ratio of the counts recorded during the main observation of 
GJ436. This was done separately for each visit. We corrected 
the counts recorded by each of the low sensitivity PMTs in the 
two pairs for all observations in a single visit using this deter- 
mined ratio. We then took the mean of the 40 Hz samples from 
the two PMTs in each pair. This procedure gave two high fre- 
quency time series, one for the "X" axis and one for the "Y" 
axis, for each visit. The following steps were done for each of 
these two sets independently. 

We attempted to correct for the inter- and intra-orbit vari- 
ances in the photometry using the measurements for the refer- 
ence stars. We first fitted a time dependent linear change in in- 
strument response to the two observations of each reference star. 
The direction of change was always the same for each reference 
star within a single visit, but not the same for all the visits. The 
slopes for the two reference stars were averaged and the result- 
ing linear trend divided from the GJ 436 observations. The mag- 



J. L. Bean et al.: A Hubble S pace Telescope transit light curve for GJ436b 



3 



nitude of this correction was less than 0.05% for 90% of the 
data points and always less than 0.1%. We repeated the analy- 
sis described in §4 on data reduced using the correction from 
only on one or the other reference star and with no correction. In 
all the cases the determined parameters were indistinguishable 
from those we present below. Therefore, our results are robust 
againsit uncertainties in this step. The correction to the GJ436 
photometry with the linear trend determined from the average of 
both reference stars yielded marginally smaller residuals from 
the light curve fit so it was utilized in the primary reduction. 

The next step was to normalize the GJ 436 data to the same 
relative flux scale. We were not successful at this no matter 
which combination of reference star data we used and the dif- 
ferences in the relative flux levels between visits were obvious 
by eye. We believe that the primary reason for this is due to the 
varying position of GJ 436 in the instrument's FOV. This hap- 
pened because the two gyro guiding mode for the HST prevented 
the same spacecraft roll angle to be used for all the visits. The 
visit groups 1, 2, and 3-6 were each obtained at different rolls. 
The FGS FOV must not have a flat photometric response at the 
level our data are sensitive to in addition to the well studied vari- 
ance in position aber ration (i.e. the optical field angle distortion, 
iMcArthur et"alll2002l) . 

Examination of the un-normalized data for the multi-visit set 
obtained at the same roll show consistent, but still not perfect, 
relative flux levels. This supports our hypothesis about the FOV 
variance. The lower level disagreement for this group could be a 
result of still small variances in FOV position even among visits 
carried out with the same roll or the effects of stellar activity 
dDemorv et al.ll200"7l) . We see no way to distinguish between the 
these two effects with the current data. 

The reference stars were never placed in the same position in 
the FOV as GJ 436, and indeed this would not have been possi- 
ble due to HST guiding and FGS pointing restrictions. Therefore, 
the reference star data cannot be used to correct the relative flux 
levels of GJ436. The time dependent response correction de- 
scribed above is still valid because that variance is most likely 
due to the thermal settling of the telescope itself and should be 
similar for all the targets. Our solution to the relative flux correc- 
tion problem was to introduce normalization parameters for the 
data obtained in each of the visits that were solved for during the 
light curve analysis described in §3. The relative flux levels for 
the visits determined from this analysis are given in Table Q] 

At this point in the reduction we had two time series (X 
and Y axes) for each of the six visits. We summed the two 
sets to make a single time series and analyzed the data as de- 
scribed in §3. This analysis yielded very poor results. The resid- 
uals were much larger than expected from counting statistics and 
clearly correlated (trends and jumps) on ~ lOminute timescales. 
Inspection of the data revealed that the source for most of the un- 
usual noise was the X axis data. When the data from the two axes 
were analyzed separately, we found that X axis data had residu- 
als twice as large as those from the Y axis despite nearly iden- 
tical count rates. Furthermore, the Y axis transit model residu- 
als do not exhibit obvious correlations like the X axis data. We 
don't have a definitive explanation for the lower quality of the X 
axis data. We note that since beginning science observations with 
the FGSlr in 2000 we have consistently (thousands of indepen- 
dent observations) obtained position residuals ~35% higher in X 
axis data compared to Y axis data when using the instrument for 
high-precision relative astrometry even though corrections deter- 
mine d from extensive calibration efforts are applied for this work 
(e.g. lBenedict et al.ll2007HBean et al1l2007l) . This discrepancy is 
similar in magnitude but opposite what was seen in data from 
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Fig. 1. Top Observed transit light curve for GJ 436b (points) with 
best fit from our analysis (line) phased to the central transit time. 
The bar on the right indicates the average error size. Bottom 
Residuals from the fit (points), which have a rms of 5.2 x 10~ 4 . 
The different point styles in both panels indicate data from the 
different visits. 



FGS3 when it was used for science observations. The effect we 
see in the GJ 436 photometry is likely related to this issue, but 
relatively larger possibly due to a lack of any sort of known ap- 
plicable correction. We ultimately decided to set aside the X axis 
data because the expected V2 reduction in counting noise from 
including these data is more than negatively compensated by the 
larger errors introduced by using it. 

The final time series that we analyzed as described in §3 was 
created by binning only the Y axis 40 Hz measurements to 60 s 
samples, which yielded 180 data points. The adopted values for 
each bin were the average of the counts and the initial error esti- 
mate was the error in the mean. These data are given in Table 2, 
which is only available electronically from the CDS. 



3. Analysis 

3. 1 . Parameter determination 

We modeled the obtained photometric ti me series of GJ 436 us- 
ing the exact analytic formulae given bv lMandel&Agoll(l2002h 
for a planetary transit. To account for the stellar limb darken- 
ing we calculated flux-weighted theoretical spectra for 18 dif- 
ferent angles from the central line of sight and integrated them 
over the unique bandpass of the FGS with the F583W filter. 
We u sed the latest version o f the PHOENIX model atmosphere 
code dHauschildt et alj 19991) for the s e calc ulations with the stel- 
lar parameters given by iBean et al.l (120061) as determined from 
a spectral synthesis analysis (T e ff - 3480 K, log g = 4.92, and 
[M/H] = -0.33). We fitted the calculated specific intensities with 
the non-linea r equation for stellar limb darkening proposed by 
IClaretl d2000l Eq. 6) to o btain the necessary lim b darkening co- 
efficients for use with the lMandel & Agol d2002) formulae {c\ = 
1.47, c 2 = -1.10, c 3 = 1.09, c 4 = -0.42). 

We used a M arkov Chain Monte Carlo (MCMC) method 
similar to that of Holman et al. (2006) to identify the best fit 
model to the data. In the analysis we fixed the values of the 
star's mass, the planet's (significantly non-zero) orbital eccen- 



4 



J. L. Bean et al.: A Hubble S pace Telescope transit light curve for GJ436b 



Table 3. Parameters for GJ 436 and transiting planet. 



Table 4. Transit times and residuals from the mean ephemeris 
for GJ 436b. 



Parameter 



Value 



Determined directly in this paper 

R* (Rq) 

Rp (R<b) 

P(d) 
in 

T c (HJD) 



505+ 0029 
4 90 +0 45 

2.643882!™°°°" 
85.80i°£ 

2454455.279241^;™°! 



Inferred 
M„ (M e ) 
P (g cm" 3 ) 



22.8+JI 



1.1 



+0.2 
-0.3 



tricity and longitude of periastron, and the radial velocity semi- 
amplitude of t he sta r due to the planet to those determined by 
Man ess et al.l d2007l) . The parameters that we allowed to vary 
were the star's radius 7?*, the transiting planet's radius R p , or- 
bital inclination i, orbital period P, central transit time at a mean 
epoch of the observations T c , and relative flux normalizations for 
each of the six visits. There were 1 1 free parameters in total. 

We adopted the standard x 2 parameter as the fit quality met- 
ric. We initially used the uncertainties in the photometry as de- 
scribed in §2, but increased them by 18% for the final analysis 
so that the best fit model had xl — 1 -0. We combined the re- 
sults from five Markov chains that were initialized with different 
randomly selected parameter values. The chains each contained 
500000 points and the first 20% were discarded to reduce the 
impact of the initial parameters. The perturbat ion sizes were ad- 
justed to give a 25% jump acceptance rate. The lGelman & Rubinl 

(fl992h Vfl statistic for the chains was within 1 % of unity for all 
the parameters, which indicates that the chain lengths were likely 
sufficient for convergence. We adopt the medians of the MCMC 
distributions as the determined parameter values. The uncertain- 
ties were estimated by the range of values that encompassed 68% 
of the parameter distributions on each side of the corresponding 
median. The model for the adopted system parameters is shown 
with the observations phased to the central transit time in Fig.Q] 
The time series has a rms of 5.2 x 10~ 4 in relative flux units, 
which corresponds to 0.6 mmag. 

To account for possible systematic errors in the determined 
parameters due to the uncertainties in the fixed parameters we 
repeated the MCMC analysis with each of these parameters 
adjusted to Icr higher and lower than their nominal values. 
We adopt the unce rtainties in the fixed parameters given by 
Man ess et al.l (120071) . We found that only the uncertainty in the 
adopted stellar mass (0.04 Mq) was a significant source of error 
in any of the parameters we determined. We added the parameter 
deviations from fixing the stellar mass to high and low values in 
quadrature with the uncertainties found from the MCMC anal- 
ysis to give the total uncertainties. The determined parameter 
values and the inferred planet mass (M p ) and density (p) along 
with their corresponding uncertainties are given in Table[3] 

The central transit time determined from the above analysis 
and given in Table[3]should be thought of as an epoch mean value 
because the observations to obtain the full light curve spanned 1 1 
of GJ436b's orbits. We carried out an analysis of the visit data 
sets separately to determine the corresponding individual transit 



T c 




O - C 


(HJD) 




(d) 


2454439.4160 ± 


0.0016 


0.0001 


2454444.7036 ± 


0.0018 


-0.0001 


2454447.3477 ± 


0.0017 


0.0001 


2454463.2109 ± 


0.0015 


-0.0000 


2454468.4986 ± 


0.0018 


-0.0001 



times. One of the visits occurred completely out of transit and 
was not included in this step. We fit a transit light curve model 
to the five visit data sets that are sensitive to the transit time us- 
ing a Levenberg-Marquardt algorithm. We allowed only the time 
for each observed transit to vary and held the other parameters 
fixed to either their previously determined or adopted values. To 
estimate the uncertainties in the transit times we used a boot- 
strap Monte Carlo method. We generated 10000 realizations of 
each visit set by randomly drawing 30 data points with replace- 
ment. We also randomly modified one of the fixed parameters 
during each iteration by drawing from a normal distribution of 
values with a width set to the parameter's uncertainty. We mod- 
ified only one of the fixed parameters for each iteration because 
the estimated errors in the parameters are correlated. We adopted 
the standard deviation in the distribution of times found from fit- 
ting the bootstrap trials as the uncertainty in the transit times. 
The determined times for each observed transit along with their 
uncertainties and deviation from the mean ephemeris are given 
in Table g] 

3.2. Assessment of possible systematic error 

The necessity of fitting for the relative flux normalizations and 
phase folding partial transit observations in our analysis is un- 
usual among the study of high-precision transit light curves. This 
raises the question of whether our results could be systematically 
wrong due to the inappropriateness of our approach. We inves- 
tigated this possibility using a combination of re-analyses of the 
data with different assumptions and simulations. 

To asses the validity of determining the physical parameters 
of GJ 436 and its planet along with the flux normalizations we 
used a simulation technique. We generated ten simulated data 
sets with model parameters from the averaged results of G07, 
D07, and S08, random noise properties similar to our observed 
data, and random offsets in the visit subsets. We fitted these 
data using the MCMC procedure described above and exam- 
ined whether there was a systematic deviation in the determined 
parameters compared to the known underlying parameters. We 
found no systematic trend in the deviations, and their distribu- 
tion followed that expected from the adopted parameter uncer- 
tainties. 

We also carried out a number of analyses on different real- 
izations of the primary data set to ascertain whether some par- 
ticular treatment of the data influenced our results. We repeated 
our analysis on subsets created by removing one of the visit sets 
to test the robustness of the results. We repeated this analysis 
six times, each time removing a different visit set. To ascertain 
whether the reduced time resolution from the 60 s binning had 
an affect, we analyzed data sets where the 40 Hz samples were 
binned to smaller ranges (10, 15, and 30 s). We also repeated our 
analysis assuming a quadratic form of limb darkening (equiva- 
lent constants to the values for the non-linear model given above 
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are j\ = 0.30, 72 = 0.49). Additionally, analyses were carried out 
by two of us separately and with slightly different methods. In 
all the cases the resulting parameters did not deviate from their 
nominal value given in Table[3]by more than 1.5cr. 

Assuming a constant orbital period, transit timing deviations 
over the course of the observations, which spanned 1 1 orbits of 
the planet, could possibly adversely affect our analysis by dis- 
torting the morphology of the folded light curve. The transit tim- 
ings we determined from our principle analysis are remarkably 
regular, but might not represent the true values as they depend on 
the physical parameters (i.e. stellar radius, and planet radius and 
inclination) determined from analyzing the phase folded light 
curve. To study this issue we repeated the transit timing analy- 
sis with the system parameters fixed to a representative of the 
G07, D07, and S08 determined values and allowing the relative 
flux normalizations to vary. From this we found possible tim- 
ing variations of up to 79 s and the x 2 for the fit was 170.9. For 
comparison the^- 2 for our principle analysis was 168.6. for 169 
degrees of freedom. 

The low value of the fit quality metric from this alterna- 
tive analysis suggests that the results from our primary analy- 
sis might not represent a unique model of the data. However, 
could true transit timings of this magnitude bias the system pa- 
rameters determined by analyzing a phase folded light curve to 
the level our results are deviant from the Spitzer results? To an- 
swer this question we generated simulated data similar to that 
described above for investigating the affect of variable flux nor- 
malizations. In this case we used transit times for the individual 
samples equal to those found from holding the system param- 
eters fixed to those determined from the Spitzer data. We then 
analyzed the simulated data in the same manner as we analyzed 
the real HST data. From the analysis of 10 simulated data sets 
we found that transit timing offsets up to at least 80 s did not 
significantly bias the results. For all the simulations the analysis 
algorithm returned the parameters that were used to create the 
simulated data within the one sigma uncertainties and without a 
systematic offset. 

The above analysis suggests that transit timings have not 
caused our results to be biased. However, this analysis still poses 
the question of whether the inverse of the previous question is 
true. That is, can adopting systematically different physical pa- 
rameters lead to spurious transit timing offsets? To study this we 
again turned to a simulation analysis. We generated simulated 
data using our determined system parameters and no transit tim- 
ing offsets. We then repeated the analysis of fixing the parame- 
ters to those determined from the Spitzer data while allowing the 
individual transit times and flux normalizations to vary. In this 
analysis we consistently found transit timing variations with the 
same magnitude (~ 80 s) as those from above, and a fit^- 2 close 
to, but still higher, than that obtained from using the correct (in- 
put) model parameters. As we know the transit timings are not 
correct because the simulated data were created assuming a con- 
stant period, we can be sure that transit timings determined from 
an analysis with assumed improper physical parameters would 
give spurious results. 

Another possible source of uncertainty in our analysis i s 
that arising from corellated, or "red," noise dPont et al.l 12006). 
We assessed the influence of this type of noise on our results 
by re-fitting different realizations of the original dat a adjusted 
with the "prayer bead" method (Mou touet al.l 12004). We used 
a Levenberg-Marquardt algorithm to perform the parameter op- 
timization. The individual data sets were modified by shifting 
the residuals from the principle analysis described in §3.1 by 
a random number and adding them back to the data. This was 



repeated 10000 times. The standard deviation of the resulting 
parameter distributions gives the uncertainties. In this investiga- 
tion we found uncertainties about half or less than those deter- 
mined from the MCMC analysis. Thus, we conclude that corel- 
lated noise is not a significant source of uncertainty in the data. 
This is likely because the data set includes measurements from 
two different orbits (i.e. the two measurements are independent 
and uncorrelated) for almost every transit phase. 

From all of the above investigations we conclude that our 
analysis is relatively unbiased, the determined parameters are a 
fair representation of the data, and the assigned parameter errors 
are realistic. 

4. Discussion 

4.1. The planet's characteristics 

Our determined mass and radius for GJ 436b is placed in con- 
text with the previous results of G07, D07, and S08, the solar 
sy stem ice giant s , and the t heoretical mass-radiu s relationships 
of lFortnevetail d2007l) and iBaraffe et al.l(l2008l) in Fig. |2] The 
radius we find differs by 1.82<x (93% significance), 1.54<x (87% 
significance), and 2.39cr (98% significance) from the values of 
G07, D07, and S08 respectively. To collapse the previous results 
into a single value we first average the G07 and S08 results be- 
cause they are from analyses of the same reduced data. Then 
we average this value with the D07 value to give R p = 4.22 R®. 
We conservatively adopt a symmetric error in this value of 0.21 
R @ , which is the largest error bar among the three results. The 
standard deviation of the three values is 0.15/? e . The radius we 
find for GJ436b is larger by 1.74<x (92% significance) than this 
representative of the previously determined values. 

The Spitzer data have an advantage in that the effect of limb 
darkening is almost non-existent due to the much longer wave- 
length of the bandpass used. However, we have carefully consid- 
ered the effect of limb darkening in our analysis through specific 
calculation of it in the FGS bandpass using a realistic model at- 
mosphere. Changing the parameters of the model atmosphere 
we used within reasonable limits or using a different functional 
form to fit the specific intensities (i.e. quadratic) does not signif- 
icantly affect our results. Therefore, our determined parameters 
are relatively robust against uncertainties in the limb darkening. 

Our data have a similar level of precision and time sampling 
as the Spitzer data. The larger uncertainties in the parameters 
from our analysis stems from the additional degrees of free- 
dom introduced by fitting for the relative flux normalizations. 
However, as we have shown above, this does not mean the re- 
sults are biased and is unlikely to be the reason for the planet 
radius discrepancy. 

The larger radius we determine for the planet cannot 
construed as a detection of the planet's effective radius de- 
pendency on wavelength arising from s pectral features (e.g. 
ICharbonneau et al.1 120021: iBarmanl l2007t iPont et al.l 120081) . al- 
though our bandpass does cover t he Na doublet near 600 nm 
(iFortnev et al .1120031; iBarmari 120071) . One reason for this is that 
we find a planet to star radii ratio consistent with G07, D07, and 
S08. If we fix the star's radius to the average of the values given 
by those authors we obtain a planet radius nearly exactly the 
same. Furthermore, our data are not sensitive to the likely size 
of the effect, which is probably at most a few percent in our wide 
bandpass. 

A plausible explanation for the differences between our re- 
sults a nd the previous one s is stellar activity. As mentioned 
above, Demor v ;t al.l d2007l) report that GJ 436 exhibits visual 
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Fig. 2. Determined mass and radius for GJ 436b from G07 (trian- 
gle), D07 (square), S08 (circle), and this paper (star). The values 
for Uranus and Neptune are indicated by the "U" and "N" char- 
acters respe ctively. The lines are theoretical relationships for ice 
giants from Fortne vet alj (120071) . where the solid line is for a 
planet with a 10% by mass H/He envelope and the dotted line 
for a planet without an envelope. The shaded region d elineates 
the possible range of values from Baraffe et al. (2008) models 
for an ice giant with a 10% H/He envelope and age ranging from 
l-7Gyr. 



flux variations with semiamplitude ~ 1%, and period ~ 50 d. 
These flux variations are consistent with periodic variations of 
emission in the Ca II H & K lines, which presumably are modu- 
lated by the star's rotation. 

Breaking the degeneracy between the planet and host star 
radii, as we have attempted to do, requ ires high-precision cover - 
age of the transit ingress and egress (ICharbonneau et aLl l2007). 
It is possible that during the different visits the planet obscured 
regions with different brightnesses arising from the spots on the 
star's surface. If this happened during the ingress and/or egress 
then the morphology of the light curve was altered and our model 
was inadequate, which could result in incorrect determined pa- 
rameters. It is also possible that if the planet obscured an un- 
usually brighter or dimmer region of the star during the other 
parts of the transit then the fitting would be negatively affected 
due to the need to phase fold the partial transit observations and 
determine the relative flux normalizations. 

Another possible explanation for the discrepancy is that we, 
G07, D07, and S08 have simply underestimated our uncertain- 
ties. In addition to analyzing the Spitzer light curve for GJ 436b, 
S08 presented an interesting study of light curves for other tran- 
siting planet systems. One conclusion he reached was that anal- 
yses of even the highest-precision data sets for a given system 
can yield significantly different parameters owing to unidenti- 
fied systematic errors in the data. In this context it is important 
to keep in mind that the G07, D07, and S08 results are not in- 
dependent as they utilize the same data, albeit reduced indepen- 
dently by two of the groups and analyzed separately by all three. 
If there is a systematic error inherent to the Spitzer data then 
the three groups would have obtained still similar, but systemati- 
cally incorrect parameters. Formally, the difference between our 
planet radius and the representative Spitzer value is significant. 
However, it could be that this difference is only representative of 
the true uncertainty in the analyses. 

In the absence of any hard evidence that our result is spuri- 
ous, we suggest that the radius of GJ 436b is larger than previ- 
ously thought. If we assume that a planet's radius dependence 



on the mass fraction of i ts envelope is linear in this regime, 
then extrapolation of the iFortnev et al.l d2007l) models to our 
measured radius gives a H/He envelope mass fraction of 15%. 
The weighted average of our radius value and the representa- 
tive value of the Spi tzer results is 4.48 + 0.17 R @ , for which the 
IFortnev et al.l d2007l) models sugg ests an envelope mass fraction 
of 10% ± 3%. The IBaraffe et all d2008j) models are fairly con- 
sistent with those of lFortnev et alj d2007|). although they im ply a 
slightly lower envelope mass fraction. IBaraffe etal J d2008) also 
emphasize the impact of the evolutionary cooling on the radius 
of a planet like GJ 436b and this adds to the ambiguity of spe- 
cific interpretations because the age of GJ 436 is quite uncertain. 
More theoretical and observational work is obviously needed to 
improve our understanding of GJ 436b and, by extension, gen- 
eral models of planet formation and evolution. 

From an observational standpoint, obtaining and analyzing 
an additional high-precision transit light curve of GJ 436b with 
a different instrument than already used would still be valuable. 
Such data would preferably be obtained in the near-infrared to 
reduce the impact of stellar variability and uncertainties in the 
limb darkening. Another interesting observational study would 
be to make a direct radius measurement of the host star with 
interferometry, similar to wh at has been done f or the transiting 
planet host star HD 189733 dBaines et al.ll2007l) . We have found 
a similar ratio of the planet and star radii as the previous studies, 
and so an additional constraint on the stellar radius would give a 
tighter constraint on the planet radius. 

What can be definitively ascertained from comparison of the 
observational results and theoretical models is that GJ 436b has a 
significant H/He envelope similar to Uranus and Neptune. Such 
similarities in structure suggest that they formed in a similar en- 
vironment. Additionally, GJ436b and the solar system ice gi- 
ants are similar in that they likely could not have formed in their 
current locations. U ranus and Neptune are th ought to be too far 
away from the Sun (Levison & Stewart 2001), while GJ 436b is 
likely too close to its host star. One hypothesis for the formation 
and evolution of Uranus and Neptune is that they were origi- 
nally in the same region of the pr otoplanetary disk as Jupiter 
and Saturn dThommes et al.l Il999l) . The later two planets ac- 
creted gas much more quickly and scattered the ice giants to- 
be outward towards their current locations, which limited their 
growth. Uranus and Neptune are therefore really gas giant planet 
cores that did not accrete gas quickly enough to complete for- 
mation. In contrast, the structure of GJ 436b likely does not re- 
quire a compa r ably v iolent past because it orbits a M dwarf. 
lLaughlin et al.l d2004l) have shown that formation of gas giants 
around low mass stars is severely hindered because these stars 
are expected to have correspondingly low surface density proto- 
planetary ^isks_and_fhedyjiamical timescale of orbiting bodies 
is longer. Lau ghlin et alj d2004l) predict a dearth of Jovian mass 
planets due to the resulting slow accretion of gas, but plenty 
of Neptune and lower mass planets around M dwarfs. In this 
context GJ 436b could be considered a failed gas giant akin to 
Uranus and Neptune. 

Furthermore, GJ436b could have experienced scattering 
due to another body like Uranus and Neptune even though 
this probably wasn't required to limit its growth into a larger 
body. The evidence for this comes from its current loca- 
tion, which is far away from its likely formation site; and 
its observed orbital eccentricity (e = 0.16), which is in di- 
rect contradiction to pre dictions of tidal circularization theo- 
ries dManess et alj |2007). A major dynamical interaction event 
after formation leading to inward scattering is an explana- 
tion that unifies both of these properties into a single evolu- 
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tionary picture. This interaction could have occurred with the 
outer object in the system that is indicated by the long term 
trend in GJ436's radial velocities. An alternative explanation 
for GJ 436b's orbital properties is disk in teraction leading to mi- 
gration and eccentricity excita tion (e.g. iGoldreich & Tremaind 
1980; IGoldreich & Saril 120031) . Continued observations of the 
GJ 436 system are needed to detect and characterize all of the 
objects it contains. Complementary theoretical studies based on 
the system census and the likelihood that GJ 436 had a low-mass 
disk are also needed to specifically assess the plausibility of the 
conceivable evolutionary scenarios. 

4.2. Transit timings 

The results of our analysis also allow us to search for deviations 
in transit times that might arise from the gravitational pertur- 
bations of another planet in the system. The times determined 
for each of the five HST partial transit observations, which are 
given in Table [4] exhibit a maximum deviation from the mean 
ephemeris of 10 s. The average uncertainty on these times is 
145 s so the deviations are fully consistent with regular transit 
times over this observational period. 

To search for possible long-term variations we compare our 
mean transit time with the published transit times from two other 
epochs. The previ ous results come fr om the initial discovery 
data presented by Gillo n et aD ((2007b), and the analyses of the 
Spitzer light curve by G07, D07, and S08 that were discussed 
above. It should be noted that the time stamps of the Spitzer 
photometry presented by G07, which was also analyzed by S08, 
are too late by 33 s (M. Gillon private communication). The 
corrected central transit times for the G07 and S08 results are 
2454280.78148 and 2454280.78174HJD respectively. We uti- 
lized these corrected times for all the invesigations described 
below. 

The transit times determined by D07 and G07 are quite con- 
sistent, but disagree with the S08 determined value by more than 
expected from the given uncertainties (~ 2cr). To compare the 
Spitzer transit times with the other data we collapse the G07, 
D07, and S08 times into a single value by first averaging the 
G07 and S08 results, and then averaging this with the D07 result. 
We adopted the largest given uncertainty from the three analyses 
(0.00016 d from D07) as the error in the collapsed transit times. 

A plot of the observed transit time residuals from our 
ephemeris (transit time and orbital period) for the three avail- 
able epochs is shown in Fig. [3] Our determined orbital period 
is consistent with, but 3 0% more precise than, the value found 
by iManess et alj d2007l) from analyzing the radial velocities of 
GJ436 only. We find a difference of 128 s between the repre- 
sentative value from the Spitzer observation and our own. The 
deviations show a nearly even trend with orbital epoch when 
also c onsidering the lower precision time given by iGillon et~aT] 
(2007b). This could be a result of an error in the period because 
the deviations are well within the range expected by the uncer- 
tainty in this parameter. Alternatively, the transit timing devia- 
tions could be due to perturbations of the host star arising from 
an additional object in t he system. 

IManess et al.l d2007l) found strong evidence for a long-term 
slope in the GJ436 radial velocities (1.4ms _1 yr _1 ), which sug- 
gests it has an additional outer companion. The nature of this 
companion is currently unconstrained, but the observed long- 
term acceleration of the star means that it has a second degree 
of motion on top of its regular orbital motion due to the known 
transiting planet. Therefore, it is reasonable, and even expected, 
that the observed transit times of GJ 436b are varying due to the 
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Fig. 3. Observ ed transit time deviat ions from our ephemeris. The 
data are from Gil lon etall d2007bi square), an average of the 
G07, D07, and S08 values (triangle), and this paper (star). For 
the later two values the error bars are smaller than the size of 
the points. The dashed lines give the ±lcr range allowed by the 
uncertainty in the orbital period propagated backwards in time. 
The dotted line shows the trend expected due to the revised value 
of the period based on fitting the transit times. 



changing position of the host star. What would be seen in such 
a case is a curvature in the transit timing offsets rather than the 
simple linear trend that would be due to an error in the orbital 
period. 

Another possible way to detect perturbations with the tran- 
sit timings is to look for deviations between the observed transit 
times and those predicted by orbital parameters found from mod- 
eling the radial velocities only. The key in this case is that the 
radial velocity orbit is relative to the system barycenter, which 
is invariant, and the transit times are relative to the position of 
the star, which should be changing. It is for this reason that we 
chose not to model our transit light curve simultaneously with 
the radial velocities. 

Rib as et all d200 8) have recently suggested that GJ 436 hosts 
an additional "Super-Earth" planet in a near 2:1 resonance with 
the transiting planet. Such a planet would likely have a strong 
influence on the observed eclipse (both transit and secondary) 
times for planet b due to dynamical interactions. In this context, 
our finding a long-term trend in the transit timings could be a re- 
sult of orbital precession caused by such interactions. Ultimately, 
a detai led dynamical study of the model proposed bv lRibas et al.l 
d2008l) in comparison with the available eclipse times and radial 
velocities should be undertaken to evaluate their claim. This is 
beyond the scope of the current paper, but we plan to present 
such a study in the future. 

If we assume that the transit timing deviations are due to an 
error in the orbital period then we may calculate a new, more 
precise value by leveraging the multiple orbits that have elapsed. 
Using our transit time and the average of the values reported by 
G07, D07, and S08 we find P = 2.643904 + 0.000005 d, which 
is an order of magnitude more precise than the value we find by 
analyzing the HST transit photometry alone. Adopti ng this pe- 
riod al so brings the earlier transit time reported by IGillon et al.l 
(2007b) in line with the other two, but this is not very significant 
due to the relatively low-precision of the value. After submitting 
this paper we becam e aware of recent transit times measured by 
lAlonso et al.l (|2008) and amateur astronomerfl These data sup- 
port our finding that the transit times for GJ436b are drifting 
from the Spitzer ephemeris. The transit times we have deter- 
mined will provide an additional benchmark for future studies 



1 See http://brucegary.net/AXA/GJ436/gj436.htm 
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aimed at detecting or setting limits on additional objects in the 
GJ 436 system. 
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